A new class of (2 + l)-d topological superconductor with Z 8 topological classification 
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The classification of topological states of matter depends on spatial dimension and symmetry class. For non- 
interacting topological insulators and superconductors the topological classification is obtained systematically 
and nontrivial topological insulators are classified by either integer or Zi. The classification of interacting 
topological states of matter is much more complicated and only special cases are understood. In this paper we 
study a new class of topological superconductors in (2 + 1) dimensions which has time-reversal symmetry and a 
7hi spin conservation symmetry. We demonstrate that the superconductors in this class is classified by Zg when 
electron interaction is considered, while the classification is Z without interaction. 



Topological states of matter (TSM) are ground states of 
gapped quantum systems which cannot be adiabatically tuned 
to a topologically trivial state without going across a phase 
transition. The first TSM discovered in nature are integer and 
fractional quantum Hall states 01 01 ■ Many topological states 
have gapless surface states on the boundary which are topo- 
logically protected within the given symmetry class. For ex- 
ample, the integer quantum Hall states have chiral fermion 
edge states which propagates unidirectionally along the edge 
of the two-dimensional system, and is robust against arbitrary 
perturbation as long as the bulk remains gapped. The number 
of chiral channels is an integer which determines the Hall con- 
ductance and characterizes the corresponding quantum Hall 
state. 

Since the recent discovery of time-reversal invariant topo- 
logical insulators (TI)(3r|5t] topological states of matter has 
been understood much more systematically in generic spatial 
dimensions and symmetries. 16rl8[] For gapped non-interacting 
fermion systems, i.e. band insulators without electron cor- 
relation and gapped superconductors in the mean-field the- 
ory sense, there are ten symmetry classes in each spatial 
dimension [loll if we only consider symmetries which pre- 
serve the positions of electrons, such as time-reversal symme- 
try and particle-hole symmetry. In each spatial dimension, 
TSM exists in five of the ten symmetry classes, which al- 
ways include three classes with integer classification, and two 
classes with Z2 classification. 0,(3] The topological classifica- 
tion depends on the spatial dimension and symmetry class in 
a periodic way[8[] which originates from the Bott periodicity. 
When electron interaction is considered, the problem is much 
more complicated and the topological classification in generic 
dimensions and symmetry classes has not been completely 
understood. Topological invariants based on single particle 
Green's function has been proposed for interacting TSMlfTll- 
[l4ll . In one dimension (ID), the classification of gapped states 
has been studied systematically by making use of matrix prod- 
uct states and entanglement spectmrn ilaJlal . In particular, L. 
Fidkowski and A. Kitaev proposed an explicit example of one- 
dimensional time-reversal invariant superconductors lfl7l for 
which the topological classification is Z for non-interacting 
states but is reduced to Zg when interaction is considered. In 
dimensions higher than one, a systematical approach of con- 
structing interacting symmetry-protected topological states 



(defined as the TSM in a certain symmetry class which be- 
comes trivial when the symmetry is broken) in spin models 
has been proposed based on group cohomology ll8ll . The re- 
lation between such states and the non-interacting TSM's in 
generic cases remains an open question. (30] 

In this paper, we propose a new class of TSM (labeled 
by D') which in two dimensions is classified by Z§. This 
symmetry class is defined by two discrete symmetries, the 
time-reversal symmetry and an additional global Z2 symme- 
try. Physically, this symmetry class consists of time-reversal 
invariant superconductors of spinful fermions with the num- 
ber of spin-up electrons (and also that of spin-down electrons) 
conserved modular 2. Without interaction, the topological 
classification is Z, the same as time-reversal breaking super- 
conductors (class D) in 2D|19]. With electron interactions, 
we show that the classification is reduced to Z§, by analyz- 
ing the edge states and topological defects. This state can be 
considered as a higher dimensional generalization of the 
topological superconductor (TSC) proposed in Ref. Il7ll . 

Definition of the symmetry class and the classification of non- 
interaction systems. Consider a 2D electron system with the 
generic Hamiltonian H = H(ci a , cJ CT ), with a =[, 4. two spin 
components, and i labels the lattice sites. The symmetries we 
consider are time-reversal symmetry defined by 



T 1 c lt T = ai, 



-Ci t 



(1) 



and the Z2 spin conservation symmetry named by R: 

R~ 1 c i ^R = Cif, R~ 1 c l iR = -en (2) 

It should be noticed that T is an anti-unitary symmetry and R 
is a unitary symmetry. In an explicit form R = (—1)^ with 
Ni the total number of spin down electrons. An equivalent 
definition of T and R symmetries which is basis independent 



T 2 = F, F 2 = 1, R 2 = 1, T~ l RT = FR 



t~Ht 



-i, R _1 iR 



(3) 



in which F = [—1) N 1-+ N ^ is the fermion number parity. The 
second line means that T is anti-unitary and R is unitary. We 
label the symmetry class defined by symmetries T and R as 
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class D' due to its relation to the class D of time-reversal sym- 
metry breaking superconductors, which will be come clear in 
later part of the draft. 

Now we consider the simple situation when the Hamilto- 
nian is quadratic in c^, cj^. Apparently, any term mixing spin 
up and down electrons such as c^Cji and c^Cji is odd un- 
der the Z2 symmetry R, such that the quadratic Hamiltonian 
#quad is decoupled into a direct sum of that of the spin up and 
down subsystems: 

#quad = #t + #4. 

H * = E ( c l> c -) h * ( 7 ) - =t < + (4) 



The single-particle Hamiltonians h\- and /i^- are related by 
the time-reversal invariance condition: 



t 



,4 



(5) 



Thus all non-interacting states in this class are direct prod- 
ucts of spin up and down superconductor ground states, and 
the spin down state is determined by the spin up state due to 
time-reversal symmetry. Consequently the topological classi- 
fication of -ffquad reduces to that of the spin up subsystem Hf, 
which is a Hamiltonian of a generic 2D superconductor, with- 
out time-reversal symmetry. (The time-reversal symmetry is 
recovered only when the spin-down state is also included.) 

It is well-known that the 2D time-reversal breaking su- 
perconductors are classified by integer, in the same way as 
integer quantum Hall states. 1 1911 The topological invariant 
characterizing this class of topological superconductors is 
the Chern number of the single-particle Hamiltonian. For 
a translation invariant system, if we denote |n,k) as the 
eigenstates of the single particle Hamiltonian h^, with cor- 
respondingly eigenvalues E n (k), the Chern number is de- 
fined as the total flux of the Berry phase gauge field in the 
Brillouin zone: C\ = J= J d 2 k(d x a y — d y a x ) with Oj = 
—i J^e <o ( n > k-l d/dki \n, k). For a system with Chern num- 
ber C\ = N, on an open boundary there are N channels 
of chiral Majorana edge states, described by the following 
Hamiltonian: 



-"edge 



\N\ 



dkvk v\ a 



(6) 



with k the momentum along the edge, and r\\ a k the edge 
state quasiparticle operator which is a superposition of elec- 
tron and hole operators, rftak satisfies the Majorana condition 
rjl ak = ??to,-fe so that the edge state fermion is its own anti- 
particle. In particular, r)^ a> k=o = T]\ a fe=o i s a Majorana zero 
mode, carrying half of the degree of freedom of a complex 
fermion ! 19l I2QTI . For N > (N < 0) v is positive (nega- 
tive) respectively, corresponding to left and right moving edge 
states. For simplicity we have set the velocity of the Majorana 
edge states to be the same. Generically the velocity can be 



different for different a but it does not change the topological 
properties. According to Eq. (0 determined by time-reversal 
symmetry, it can be shown straightforwardly that and h/ 
have opposite Chern number. Correspondingly the edge states 
of them consist of the same number of Majorana edge states 
with opposite chirality. For C\ = — C\ = N > 0, the edge 
theory of the whole system is 



Hedge — 



\N\ 

dkvk [v\akmak 
a=l 



VlakViak 



(7) 



Since Tj-aw is a superposition of c^m, the action of R sym- 
metry on T^f(4.) fc is R^ 1 rj tak R = r]^ ak , i? _1 ^ a fci? = -174.0k- 
Consequently, no mixing term such as vlak^lak' is allowed in 
the Hamiltonian H c a gc if we consider generic perturbations 
preserving the R symmetry. 

In summary, from the analysis above we conclude that the 
2D superconductors with T, R symmetries are classified by 
integer, and the class labeled by integer N has an edge the- 
ory of \N\ flavors of non-chiral Majorana fermions, which are 
protected to be gapless by the two symmetries. 

Classification of interacting theory. Now we consider 
the effect of electron interaction in this system. With in- 
teraction the terms mixing spin up and down electrons are 
allowed by the R symmetry, such as a pair hopping term 

H p h = J2ijkt (tijkicl^ctf-ckicii + h.c)j with h.c. denotes 
the Hermitian conjugate of the first term. The time-reversal 
symmetry leads to some requirements on the matrix elements 
UjH = tkiij- Other terms such as cj^c]^c^c^ can also be 
considered. To see the consequence of such interaction terms 
on the Z topological invariant, it is most convenient to study 
the stability of the edge states described by the Hamiltonian 
<J2J under interaction. In the following we will show that the 
Hamiltonian (|7]) with N = 8 can be gapped by an interact- 
ing term without breaking the T, R symmetries, so that the 
classification of the bulk topological states can at most be Zs. 

This conclusion is obtained based on the results of Ref . lfl7ll 
in (l + l)-dimensions. It was shown explicitly in Ref. Il7ll that 
the following interacting term induces a mass to 8 channels of 
free Majorana fermions given by Eq. (0: 



Hint = 



Habcd — 



Hi234 + -^5678 + -"1256 + #3478 + #3456 
— #2367 — #1458 — #2358 — #1467 
+#1357 + #2468 — #1368 — #2457 

dxMe abcde ^h Vet71ftlU7lH 



#i 



278 



(8) 



For 8 Majorana fermions with the same velocity, there is an 
SO (8) symmetry between the Majorana fermions generated 

by the operators J ab = j- ([r/ot,%|-] + bz<4, TfclJ)' and the 
Majorana fermion operators i] a form vector representation of 
this 50(8). The interaction term given above breaks the 
50(8) symmetry to 50(7), but the 50(7) subgroup is not 
the one which preserves some vector n a of 50(8), but the 
one which preserves some spinor ip a of 50(8). In the supple- 
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mentary material we provide an explicit explanation why the 
interaction term (0 leads to such a nonconventional symmetry 
breaking. 

It is straightforward to verify that H int defined in Eq. © 
preserves both T and R symmetries. Thus the existence of 
such a mass term suggests that the quadratic Hamiltonian © 
with Chern number N = 8 for spin up electrons is topologi- 
cally trivial, since its edge states are not topologically robust 
when interaction is considered. For N > 8 one can add the 
same mass term for any 8 of the Majorana fermions, leading 
to the conclusion that the Hamiltonians with Chern number 
N and N — 8 have topologically equivalent edge states. If the 
edge states for N < 8 is stable, one is lead to the conclusion 
that the topological classification of T, R invariant topological 
superconductors is Z 8 rather than Z. 

Edge theory and the relation to Id BDI class Now we present 
further evidences that the systems with Chern number TV = 
0, 1, 7 are indeed topologically distinct from each other. 
Consider to add the following spatially inhomogeneous mass 
term to the edge states: 



H m = I dx 



\N\ 



im{x)r]^ a (x)r]i a {x) 



(9) 



with m(x) G R real. Physically, such a mass term can 
be obtained by coupling the edge states of the TSC with 
an s-wave superconductor by superconducting proximity ef- 
fect. In the lattice Hamiltonian the s-wave pairing term 

has the form of Ha = J2i (^ c !t c ll h.c.\, which pairs 



each fermion Cj CT with its time-reversal partner T 



rT = 



aci^ a . Such a property of s-wave pairing is generic and ba- 
sis independent ly ill , so that in the basis of low quasiparticles 
the pairing term has the form of ([9). The mass term m(x) is 
proportional to the imaginary part ImA of the s-wave pairing, 
since the real part ReA preserves time-reversal symmetry and 
does not induce a mass of the edge states. 

In particular, consider the configuration of m(x) with a do- 
main wall where m{x) changes its sign, as shown in Fig. Q] 
(a). This configuration corresponds to a Josephson junction 
between two s-wave superconductors in proximity with the 
TSC, as shown in Fig. |T](a). As is well-known, such a mass 
term domain walHeads to Majorana zero modes localized on 
the domain wall|22[|. For Chern number N there are |JV| zero 
modes on the domain wall. As is known for other symmetry 
classes, such as the quantum spin Hall state with T-breaking 
domain wall on the edge |23Q , such zero modes on the domain 
wall of a symmetry breaking mass term can carry fractional 
quantum numbers, in which case it can be used as a probe 
of the topological state even if electron interaction is consid- 
ered. For example, in the case of quantum spin Hall insulator 
the T-breaking mass domain wall on the edge traps a frac- 
tional chargeJZj. In the current system, to see if the zero 
modes carry fractional quantum numbers, one can first ana- 
lyze the symmetry of the system with the domain wall. The 
mass term H m in Eq. (O breaks both T and R symmetries, 



but preserves the combined symmetry T = T o R. From the 
definition of T and R in Eqs. (Q~|i and (f2j) one can see that 



"it 



T 



"a 



T 



Cif; 



so that T 2 = 1. Thus T 



can be considered as the time-reversal symmetry for spinless 
fermions, which is the symmetry defining the BDI symmetry 
classUQl. 

The domain wall zero modes are thus topolog- 
ically equivalent to the edge zero modes of a Id topological 
superconductor in the BDI class lfl7ll . The relation between 
the domain wall zero modes and the edge state of Id TSC 
can be seen most explicitly in the geometry shown in Fig. [T] 
(b). A strip of 2d TSC with the upper and lower surface states 
gapped by opposite mass m and — m can be considered as a 
Id TSC in the BDI class. With open boundary in the hor- 
izontal direction, the end of the Id system becomes a mass 
domain wall between m and — m which has Majorana zero 
modes if the topological invariant is N ^ mod 8. In Ref. 
I BE 

12411 it was shown that the edge zero modes of a ID 
system are classified by projective representations of the sym- 
metry group, or equivalently, the second cohomology of the 
symmetry group. Here the symmetries of the system are T 
and the fermion number parity F which is always a symmetry 
for fermion systems. T and F commutes with each other, and 
form the symmetry group of Z2 x Z2. There are 8 distinct 
projective representations of this group which corresponds to 
N = 0, 1, 7 number of Majorana zero modes. Therefore 
from the results on ID BDI class one can conclude that the 
domain wall zero modes for N = 1, 2, .., 7 is stable, so that 
the corresponding bulk states are topologically distinct. 

It is helpful to make more comments on the relation of the 
N = 8 edge state problem with the (1 + l)-d topological su- 
perconductor studied in Ref. Il7ll . We start by considering 
more general mass terms in the edge theory which preserves 
the T symmetry. A quadratic mass term is allowed in the edge 
theory in the form of H m = if dxA a ii] a ^(x)r]f ) i(x). The re- 
quirement of Hermitivity of the Hamiltonian and T symmetry 
requires the condition A a i, = Af, a € R, i.e. the matrix A is 
real and symmetric. Such a massive Majorana fermion can be 
viewed as a (1 + l)-d topological superconductor in BDI class 
studied in Ref. 11711 . with a topological invariant N = Ind(^4) 
defined as the number of negative eigenvalues of mass matrix 
A. If only quadratic mass terms H m are considered, the topo- 
logical classification would be Z since some eigenvalues of 
A a i, must vanish at the transition between states with different 
index N, leading to the topological phase transition. However 
it is shown in Ref. Jl7ll that a quartic mass term H- lnt given in 
Eq. ^ is possible for N = 8, such that it is possible to have 
a cross over between two A a f,'s with index different by 8. For 
example, consider the mass term 



H{6) 



8 

; cos 6 / dx ^2 Va^Vai + sin 6H int ( 1 0) 



For m > 0, the Hamiltonian at 8 = 0,tt has index and 8, 
respectively. By adding the interaction term _ffi n t, one obtains 
an interpolation between H(Q) and H(ir) without closing the 
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(a) m(<). 



(b) 




FIG. 1: (a) Illustration of the T-breaking mass domain wall along 
the edge of the 2d TSC. The mass is induced by a proximity effect 
of the TSC with an s-wave superconductor, and the mass domain 
wall corresponds to a Josephson junction between two such s-wave 
superconductors (see text), (b) Illustration of a strip of 2d TSC with 
the edge states on the two edges gapped by opposite mass terms m 
and — m, which is topologically equivalent to a Id TSC in the BDI 
class, with Majorana zero modes at the end. (c) Illustration that a 
sphere with two topological defects is equivalent to a cylinder with 
periodic boundary condition. 



gap. On the edge of (2 + l)-d topological superconductor, the 
additional R symmetry requires the quadratic mass term A a t 
to vanish. Thus the edge state theory can be viewed as the 
topological phase transition theory between different topolog- 
ical states in the (1 + l)-d BDI class. Due to the presence of 
the quartic mass term H- mt , there are only 8 distinct phases in 
(1 + l)-d BDI class. Thus there are 8 distinct phase transition 
theories at presence of the R symmetry, which becomes the 
robust topological edge states in the (2 + l)-d case classified 
by Z 8 . 

Topological defect in the bulk. Another important topolog- 
ical property of the topological superconductors in D' class 
is the topological defects in the bulk. Generically, supercon- 
ductors have vortices as point-like defects. If we restrict the 
discussion to T, R invariant defects, the only vortices pre- 
serving the symmetries are those with opposite vorticity for 
spin up and down components. Majorana zero modes can 
appear in the vortex core. To understand the zero modes 
in the vortex core, one can relate the vortex zero modes to 
the edge states. As shown in Fig. Q] (c), a sphere with two 
holes is topologically equivalent to a cylinder with two edges. 
The edge state on one edge of the cylinder is described by 
the edge Hamiltonian (O compactified on a finite size circle: 



N\ 



,t 



,t 



Their are 



two possible boundary conditions for the Majorana fermions 
on the cylinder, periodic or anti -periodic. When the cylinder 
is deformed to the sphere with two holes, the periodic and 
anti-periodic boundary conditions are mapped to the sphere 
with and without a vortex in each hole, respectively. For 
periodic boundary condition, the wavevector k takes values 
of k = 2irn/L, n G Z with L the perimeter to the edge. 
When the perimeter L shrinks, the states with nonzero k has 
higher and higher energy oc 1/L, and the only important low 



energy modes are the Majorana zero modes if\- a o , 774,00 with 
k = 0. The coupling term between different spins such as 
i'mabV^aoVibo is forbidden by the T and R symmetries. How- 
ever, different from the case of one-dimensional edge states, 
the zero modes of the same spin can be coupled to open a gap. 
For example if there are two zero modes r7-|-(n a o, a = 1,2, 
a mass term im (?7tio?7t20 — r liw r ll2o) can be added to com- 
pletely gap out the vortex zero modes. Therefore the vortex 
zero mode is classified by Z2 although the bulk topological 
invariant is Zg. This is similar to the vortices in the D class 
which is also classified by Z2. When there are odd pairs of 
Majorana zero modes in the vortex core, no mass term can 
be introduced, and at least one pair of Majorana zero mode 
remains stable. Under time-reversal T the single pair of Ma- 
jorana zero mode forms a Kramers doublet, in the same way 
as the vortices in Dili class with only T symmetry 12511 . 
Conclusion and discussion. In summary we have demon- 
strated that in (2 + l)-dimensions there is a class D' of TSC 
defined by two discrete symmetries which has a Zg topolog- 
ical classification. We studied the edge and bulk topological 
defects of the D' TSC. The (2 + l)-d D' TSC is related to 
the (1 + l)-d BDI class of TSC with symmetry f by adding 
one spatial direction and one additional discrete symmetry T. 
Such a relation between TSM in different dimensions applies 
to more generic cases. In Kitaev's "periodic table" in Ref. 
IH] it is shown that the topological classification remains the 
same if the spatial dimension is increased from d to d + 1 
while adding one symmetry T satisfying T 2 = — 1. Our re- 
sult demonstrates that such a relation between TSM in differ- 
ent dimensions with different symmetry can be generalized to 
interacting systems. The physical reason of such a correspon- 
dence has been illustrated in Fig. (Q~|) (b). Consider a d + 1- 
dimensional system with one of the dimensions (the vertical 
direction in Fig. [U(b)) small compared to other dimensions. 
When the added symmetry (T in the case of D' class) in the 
d + 1 dimensional system is broken on the boundary but pre- 
served in the bulk, generically the surface states are gapped 
and such a d + 1-d system can be viewed as a d-dimensional 
TSM in the known symmetry class. However there are always 
two ways to gap the surface states, which corresponds to do- 
ing a T transformation to the boundary symmetry breaking 
terms. A domain wall forms between the two regions with 
opposite symmetry breaking mass terms, as is shown in Fig. Q] 
(b). If the d-dimensional TSM obtained in this way is topolog- 
ical nontrivial, the domain wall must carry zero modes which 
are topologically equivalent to the boundary states of the d- 
dimensional TSM. 

We would like to emphasize an important difference be- 
tween the cases of interacting and non-interacting systems. 
For non-interacting TSM, adding a symmetry with T 2 = — 1 
has the same effect as removing a symmetry with T 2 = 
— 10|. In the example of Id BDI class, if one removes the 
T symmetry one obtains the D class (generic superconduc- 
tors). Instead, adding the symmetry T leads to the D' class 
we discussed. Both D and D' class TSC are classified by in- 
teger in the non-interacting classification, but they are distinct 
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in interacting systems. This is consistent with the fact that the 
edge domain wall construction only works when a discrete 
symmetry is added, rather than removed. More systematical 
discussion on this dimensional reduction approach to TSM is 
reserved for future works. 

I acknowledge helpful discussions with Lukas Fidkowski 
and Alexei Kitaev. This work is supported by the Packard 
Foundation. Recently I became aware of several indepen- 
de nt appro aches to similar 2d topological states classified by 
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Appendix: Further analysis on the interaction Hamiltonian © 



In this appendix we provide an explicit explanation on the symmetry property of the interaction term ([8]). The analysis is 



already presented in Ref. II 1711 . but we would like to provide an alternative illustration that is hopefully simpler and more 
explicit. 

To see the symmetry breaking induced by the interaction term (0, we introduce the Clifford algebra Cliff (8,0) with the 
generators r o ,a = l,2,...,8 satisfying 

{r Q ,r b } = 25 ab (A-i) 

r a can be represented by Hermitian matrices. If we don't require the representation of Clifford algebra to be real, the minimal 
dimension of faithful representation for the algebra Cliff (2n, 0) is 2". If we require the representation to be real T a = TJ, usually 
the minimal dimension is larger than 2™. The special property of Cliff (8, 0) is that the faithful representation with minimal 
dimension 2 4 = 16 is also real. To obtain an explicit understanding, one can consider the following matrix representation of T a : 







r 2 




r 3 


= Vy<J X Tyll Z 


r 4 


Z/y G x Tyflx 


r 5 


= -v y a z T x [L v 


r 6 


= Vy&zTylM) 


Tr 


= Vy(T z T z Hy 


r 8 





(A-2) 

Here k$, Oj, Tj, fa are matrices in four independent 2x2 spaces, which are defined as the Pauli matrices for i = x,y,z, and 2x2 
identity matrices for i = 0. The multiplication in the equation above shall be understood as direct product, so that T a are 16 
matrices. It can be directly verified that this representation of T a 's is real, and satisfies the Clifford algebra defined in Eq. (IA-11) . 

We define the generators T ab of vector representation of 50(8) as 

[T ab ] cd = -i (S C X ~ S d J c b ) (A-3) 

A generic 8x8 orthogonal matrix is written as O = cxp \iT ab 9 ab '\. The spinor representation of 50(8) is generated by 
the commutators T ab = [r a ,r o ] /4i. The elements of the spinor representation are U = cxp [ir a f,6' afc ] . T a carries the vector 
representation of 50(8) under the similarity transformation defined by T a — > U~ 1 F a U = O ab T b . In the choice of representation 
( 1A-2I >. all r a b's have the form of T ab = vq... or r ao = v z .... Consequently all r ao are block diagonal with two 8x8 blocks. In 
other words, the spinor representation is reducible. The two 8x8 blocks are the spinor representation and the conjugate spinor 
representation of 50(8). Moreover, because T a is real, T ao is purely imaginary and the rotation in the spinor representation U 
is real. If we denote 

T ab = ( 7 <f ° ) (A-4) 



7 a " 

one can define the mappings from vector representation to the spinor representations 

V+ : O = exp [iT ab 9 ab ] — > U+ = cxp [ij+8 ab ] , 

if- : O = exp [iT ab 9 ab ] — > U- = cxp [i% b 6 ab ] (A-5) 

Because U± are also real orthogonal matrices, the two maps <p± are isomorphisms between the vector representation and the 
spinor representations. This is the so-called triality property of 50(8). 

Now we consider 8 Majorana fermion operators rj a . The Hilbert space of such as system is 16 dimensional, the same as that 
of 4 complex fermions. Since rj a satisfies the Clifford algebra in Eq. (IA-11) . we can take a basis in the Hilbert space so that 
the matrix representation of operators rj a is T a . The time-reversal symmetry T acts on r/ a as T^ 1 rj a T = r\ a - Therefore in the 
representation given in Eq. ( IA-21 ) one can take the time-reversal symmetry to be T = K which is the complex conjugation. Now 
we want to define a time-reversal invariant Hamiltonian in this Hilbert space with a unique ground state. The simplest choice is 
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a projection operator: 



proj 



( 



\ 



(A-6) 



which has all entries vanishing except the lower right corner. Such a projector Hamiltonian breaks the 50(8) of the lower 8x8 
block to 50(7). However, this 50(7) is the 50(7) which rotates the first 7 columns of the lower 8x8 block. This is not 
the usual 50(7) subgroup of 50(8) which preserves a given vector. Instead, it's the subgroup which preserves a given spinor 
(0,0,....,!). 



We would like to express such a Hamiltonian in the Majorana fermion operators r] a , i.e., T a in this representa- 
tion. In general, a complete basis of 16 x 16 real matrices can be obtained by multiplying different r a 's. The set of 
{1, T a , r a rh, r a r;,r c , rilV.-.Ts} (with all a, b, c... in each term different from each other) contains exactly 



8 

E 

71=0 



independent real matrices. Because T a are all off-diagonal, the multiplication of even l? a 's are non-vanishing in the diagonal 
8x8 blocks and those of odd IVs are non-vanishing in the off-diagonal 8x8 blocks. In this set, the matrices with nonzero 
overlap with flp ro j are those which are both symmetric and block-diagonal. Beside the identity 1, there are only two set matrices 
satisfying these two conditions: TaT^rcr^ and F = rilV-.-IV In the representation ( IA-2t F = —v z . Physically F is the 
fermion number parity. Since the Hamiltonian ( IA-6b vanishes in the upper 8x8 block, we can focus on the basis of the lower 
8x8 block obtained by the projection: 



1 



1 



labcd = +F )^l E 



abcdef gh 



(A-7) 



The antisymmetric tensor e abcde fa h i s introduced to extract only the components of Tar^FfT^ with a, b, c, d all different from 
each other, labcd satisfies the orthogonality condition 



1 



2 Tr [labcdlefgh] 



l + F 



abcdij kl , 
^ t-ijklefgh ~> ^abcdefgh 



(A-8) 



For example 71234 = |(1 + F)r 5 r 6 r 7 r 8 = \ (rir 2 r 3 r 4 + r 5 r 6 r 7 r 8 ), which is orthogonal to other j abcd except 7 5678 = 
71234- The number of independent labcd is therefore ^f^T = 35 which is consistent with the number of symmetric traceless 



8x8 matrices 35 = — 1. Using the orthogonality condition one can expand the Hamiltonian to 



-Hproi — 



l + F 
16 



-tr 



l + F 



/ / gtr [jabcdHpro]] 1 labcd 
abed 



(A-9) 



with the first term taking care of the trace part and the second term an expansion of the traceless part in jabed- Here the sum 
Y^'abcd stands for the sum over the 35 independent groups of abed. For the particular Hamiltoinan H pro j, we have 



tr [■JabcdHproj] = — [labcd] 



16,16 



(A-10) 



which is the last diagonal component of labcd- In the representation given by Eq. ( IA-2b . labcd is also a direct product of 
Pauli matrices or identity. Consequently the 16, 16 component is only nonzero if labcd — Vi<jjTkHi with all k, I — or z. 
Therefore there are only 7 = 2 3 — 1 terms with nonzero contributions. They can be explicitly found as 



71234 = 
71357 = 



''(I 



v z v 

CTzT-o^o, 71256 = 



-0- 2 ToMz, 71368 



-CT Z T Z ^ , 71278 



-O0T ^ z , 71458 



2 

- v~ 



-cr T z fi a 



-&0T z fJ-z, 71467 



^0 



: cr z r z /i z (A-ll) 



8 



Therefore we obtain 

1 + F 1 

-ffproj = g [71234 + 71256 + 71278 + 71357 ~ 71368 ~ 71458 ~ 71467] (A- 12) 

jabcd can be written in product of r a 's, which in the operator form are rja's. For example 71234 = \ (171*72773174 + TW)§rftr\%). 
Compare Eq. (IA-12b with Eq. ((SJ one finds that except for a trivial constant term, they are identical if we only take the zero 
modes r\ a ,k=o\{X) m Eq- ®- Consequently, we have demonstrated that the Hamiltonian ((HJ leads to a unique ground state in the 
Hilbert space of the 8 Majorana zero modes. 



